Dynamics in two-leg spin ladder with a four-spin cyclic interaction 
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We study two-leg Heisenberg ladder with four-spin cyclic interaction using the (dynamical) 
density-matrix renormalization group method. We demonstrate the dependence of the low-lying 
excitations in the spin wave, staggered dimer order, and scalar-chirality order structure factors on 
the four-spin cyclic interaction. We find that the cyclic interaction enhances spin-spin correlations 
with wave vector around momentum {qx, qy) = (-1,0). Also, the presence of long-range order in the 
staggered dimer and scalar-chirality phases is confirmed by a (5-function peak contribution of the 
structure factors at energy uj = Q. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.40.Gb, 75.40.Mg 



For many years, it had been generally believed that the 
magnetic properties of undoped high-Tc materials can 
be well described by the two-dimensional (2D) Heisen- 
berg model with nearest-neighbor exchange interaction 
J. However, four-spin cyclic interaction K has been in- 
creasingly recognized as a non-ncgligiblc correction to the 
Heisenberg model. The cyclic interaction comes from the 
forth-order processes in the strong-coupling limit of the 
single-band Hubbard model at half fillingfi The impor- 
tance of this interaction was initially proposed in the 2D 
solid ^He, which has the hard-core correlations between 
spin-i fermionsi^ 

In fact, a substantial value K — 0.24J was proposed 
for 2D copper oxide La2Cu04 by an accurate fit of the 
magnon dispersioni^ A close value was also suggested by 
an analysis of the Raman-scattering datai^ Such mag- 
nitude of the four-spin cyclic interaction must have a 
considerable influence at least quantitatively on the low- 
energy spin physics. Similar situations have been re- 
ported for several two-leg spin-ladder systems: the ex- 
change interactions were estimated as J\\ = J± = 110 
meV and K = 16.5 mcV for La6Ca8Cu2404i (Ref. 0); 
J|| = 186 meV, Jj_ = 124 meV, and K = 31 meV for 
La4SrioCu2404i (Ref.[i); Jy = 165 meV, = 150 meV, 
and K = 15 meV for SrCu203 (Ref. 0), where J|| and J± 
are exchange interactions in the leg and rung directions, 
respectively. 

Motivated by those observations, the ground-state 
properties of two- leg spin-i Heisenberg ladder with the 
four-spin cycHc interaction have been intensively stud- 
iedi^'^'^ i^^i^^ Also, the effect of magnetic field on the 
ground state has been investigated ] Furthermore, the 
spectral features of the spin structure factor have been 
examined by the exact diagonalization, perturbation the- 
ory, and density-matrix renormalization group (DMRG) 
methodF^iii^ii^iii The spin dynamics for small cyclic in- 
teractions is thus well understood, but the dynamical 
properties for other correlations and/or relatively large 
cyclic interactions are still open. In this paper, we study 
the dynamical structure factors of staggered dimer order, 




FIG. 1: Lattice structure of the two-leg Heisenberg ladder. 
J|l (Jx) is the exchange interaction in the leg (rung) direction 
and K is the four-spin cyclic interaction. The x- (y-) axis is 
defined as the leg (rung) direction. 



scalar-chirality order, and spin waves for a wide range of 
the four-spin cyclic interaction to give a deeper insight 
into our knowledge of the low-lying excitations, using the 
dynamical DMRG (DDMRG) method.-i^ 

The Hamiltonian of the two- leg spin-i Heisenberg lad- 
der with the four-spin cyclic interaction is given by 

x.,y X 

+kY,{p,. + p~^), 



(1) 



with the cyclic permutation operator 



Px + Px"^ ^ ^x,l ■ Sx,2 + Sx+lS ■ >5'a;+i^2 + 5^,1 ' Sx+1,1 

+ Sx,2 ■ Sx+1,2 + Sx,l ■ Sx+1,2 + Sx,2 ■ 'S'a- + l,l 

+4(S'a:,l • Sx,2){Sx + l,l ' S'2; + l,2) 

+'l{Sx,l • Sx + l,l){Sx,2 ■ S':i- + i_2) 

— M^x,! ■ Sx + l,2)iSx,2 ■ Sx + l,l) (2) 

where Sx.y is a spin-i operator at a site {x, y) [sec Fig.[T]. 
For simplicity, we focus on the case of J|| = J± ~ J and 
take J = 1 as the unit of energy hereafter. The ground- 
state phase diagram was obtained in Ref. H as follows: 
the system has a rung singlet phase for —3.33 < K < 
0.23, a staggered dimer long-range order (LRO) phase 
for 0.23 ^ K ^ 0.5, a scalar-chirality LRO phase for 
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0.5 < K < 2.8, a dominant vector chirality phase for 
2.8 < and a ferromagnetic phase for K < —3.33. 
Let us define the dynamical structure factor as 

x5{uj-E, + Ea), (3) 

where is the v-ih. eingenstate with the eigenenergy 
En and is the Fourier transformation of the quantity- 
dependent operator Op. The (5-function is replaced by 
a Lorenzian with width rj in our numerical calculations. 
We now study the following three kinds of the dynamical 
structure factor corresponding to three phases at —3.3 
K 2.8. The first is spin structure factor S{q,u!) with 
the operator 

Or = Sly, (4) 

where y is z-component of the total spin, the second 
is dimer order structure factor D{q,Lj) with 

Or — Sx — \^y ■ Sx.y Sx.y ' Sx-\-\^y: (5) 

and the third is scalar-chirality structure factor C{q,uj) 
with 

Op = Sx,l ■ {Sx+1,1 X S'a;-|-i,2)- (6) 

By integrating Eq. wc can easily obtain the static 
structure factor 

A{q) = (^o|0-AHV'0>. (7) 

We employ the DDMRG methodic which IS an ex- 
tension of the standard DMRG methodJ^ It has been 
developed for calculating dynamical correlation func- 
tions at zero temperature in quantum lattice models. 
This method has been successfully applied to the one- 
dimensional Heisenberg modeli^ We now calculate the 
dynamical structure factor ^ with applying the peri- 
odic boundary conditions in the leg (x) direction. We 
fix the system length L ~ 32 and 77 = 0.1 if not other- 
wise stated. In the DDMRG calculation, a required CPU 
time increases rapidly with the number of the density- 
matrix eigenstates so that we would like to keep it as few 
as possible; meanwhile, the DDMRG approach is based 
on a variational principle so that we have to prepare a 
'good trial function' of the ground state with the density- 
matrix eigenstates as much as possible. Therefore, we 
keep m — 600 to obtain true ground state in the first 
ten DDMRG sweeps and keep m = 300 to calculate the 
spectral functions. In this way, the maximum truncation 
error, i.e., the discarded weight, is about 3 x 10~^, while 
the maximum error in the ground-state and low-lying ex- 
cited states energies is about 10~^. 

To begin with, we consider the spin structure fac- 
tor. In Fig. [21 we show the DMRG results of the static 
and dynamical spin structure factors for the rung sin- 
glet {K = 0.1), staggered dimer LRO {K = 0.4), and 
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FIG. 2: (a) Static spin structure factor, (b) Dynamical spin 
structure factor for K ^ 0.1 (Top), K = 0.4 (Middle), and 
K = 0.7 (Bottom). Left and right panels correspond to the 
results for qy = and qy = tt, respectively. Insets: lower- 
edge of the two-spinon continuum {qy = 0) and one-spinon 
dispersion (qy = tt). The dashed line denotes the perturbative 
result uj{qx,qy = tt) = 1.186-f 0.558 cos(ga;) —0.271 cos(2gi:) -|- 
0.071 cos(3g^). 



scalar-chirality LRO {K = 0.7) phases. The dispersion 
relations i-u{q) are also plotted in the insets. The spectra 
for qy ~ and n exhibit the two-triplon and one-triplon 
contributions, respectively. In the absence of the cyclic 
interaction, i.e., K = 0^^^ it is known that the spin dis- 
persion has two minima at qx = 0, tt and a maximum at 
qx ~ 27r/3 for qy = 0; whereas, two minima at qx = 0,7r 
and a maximum at q^, 7r/3 for qy = tt. Those features 
have been confirmed to remain qualitatively unchanged 
at if < 0.075.^ For K = 0.1, however, the minima at 
{Qx,qy) — (tt, 0), (0, tt) are no longer visible [see the in- 
sets of the top panels in Fig. [5]Jb)], i.e., the dispersions 
are nearly flat around {qx,qy) ~ (tt, 0), (0, tt). The one- 
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triplon excitation {qy = tt) is in good agreement with the 
perturbative result This is also consistent with other 
numerical study.— 

When the cyclic interaction is further increased to 
K = 0.4, we can see a drastic change of the spectra 
for both Qy ~ and tt: especially, an enhancement 
of peaks around {q^^qy) = (§,0) and a reduction of 
peaks around {qx,qy) = (tt, tt) are derived. It is be- 
cause the cyclic interaction leads to a repulsive inter- 
action between neighboring rung triplets. In addition, a 
node emerges at qx = ^ for both qy values and a re- 
lation U!{qx,qy = 0) = Lu{7r — qx,qy = tt) appears to 
be satisfied. They would indicate a two-fold degener- 
ate ground state with a broken translational symmetry, 
which is consistent with the staggered dimer order state. 
For K = 0.7, the peaks around {qx,qy) = (f ,0) are still 
more enhanced and, whereas, the low-energy spectral fea- 
tures for qy = TT seem to be much reduced. Actually, the 
one-triplon contribution is shunted off to the high-energy 
excitations since the static structure factor for qy = tt is 
not much suppressed. In short, from the standpoint of 
spin-spin correlation the four-spin cyclic interaction may 
work for enhancing a spin-density wave with wave vector 
{Qx,Qy) = (fjO) and for reducing the antiferromagnetic 
correlation S'(7r, tt) [see Fig. [21[a)]. 

Next, we turn to the dimer-order structure factor. Fig- 
ure [3] shows the DMRG results of the static and dynam- 
ical dimer order structure factors for K ~ 0.1, 0.4, and 
0.7. For comparison, the results of D{q,u!) ioi K = 
arc shown in the insets of the top panels of Fig. [Sj^b). 
In the rung-singlet phase, the ground state is approxi- 
mately expressed as the product of local rung singlets 
with gap A ^ 0{J±). The lowest excitation comes 
from the formation of a leg singlet with coupling en- 
ergy ^ as well as the collapse of two rung singlets. 
For K = 0, therefore, undispersive sharp peaks appear 
around u) = 0{2A — ^) ^ 1.5 for qy = 0; whereas, 
the spectra for qy ~ tt consist of broad continua at 
uj > 0{2A - J„). 

When small cyclic interaction (K = 0.1) is introduced, 
we can see a strong influence on the continua around 
ilxTQy) = (''r,7r), i.e., they are significantly shifted to- 
wards lower energies. It implies that the gap A is re- 
duced rapidly as K increases. For K = 0.4, the continua 
are further drastically changed: a pronounced peak ap- 
pears at {qx,q_y) = (tt, tt) and a; ~ 0; also, most of the 
spectral weight concentrates around the peak. On the 
other hand, the spectral weights for qy = Q are totally 
suppressed. The pronounced peak is well fitted by a 
Lorcnzian with rj = 0.1, as shown in the inset of the mid- 
dle panel of Fig. ^h). In other words, the spectrum for 
{Qx, Qy) = (i", tt) may consist of a (5- function peak at w = 
and a gapfuU continuum. They would be a signature of 
long-range staggered dimer order. For K = 0.7, the spec- 
tral weights around {qx,qy) = (7r,7r) are much reduced 
and a gap opens. It means that the staggered dimer 
order is no longer dominant in the ground state. Never- 
theless, the spectral weights around {qx,qy) = (t^jT^) are 
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FIG. 3: (a) Static dimer-order structure factor, (b) Dynam- 
ical dimer-order structure factors for K = 0.1 (Top), K = 0.4 
(Middle), and K = 0.7 (Bottom). Left and right panels cor- 
respond to the results for qy = Q and qy = tt, respectively. 
Insets of the top panels: D{q,uj) for K = Q with L = 16 and 
77 = 0.2. Inset of the muddle panel: A Lorenzian fit of the 
peak at {qx,qy) = {tTjTt) and uj ^ with rj = 0.1. 

still significant, as seen in Fig. ^a), so that the dimer- 
order correlation could just be changed from long-range 
order to short-range order. It is consistent with the fact 
that the staggered dimer order parameter is finite even 
in the scalar-chirality LRO phased 

Finally, we look at the scalar-chirality structure fac- 
tor. The DMRG results of the static and dynamical 
scalar-chirality structure factors ior K = 0.1, K = 0.4, 
K = 0.7 are shown in Fig.[31 For K = 0.1, the lowest ex- 
citations are described by (almost) undispersive peaks 
around u> ^ 0{2A — -j-) in analogy with the dimer- 
order structure factor. For K ~ 0.4, the spectra form 
a continuum bounded by the branches uj{qx) ~ Asm{qx) 
and uj{qx) 2Asm{qx/2), except that a gap opens at 
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FIG. 4: (a) Static scalar-chirality structure factor, (b) Dy- 
namical scalar-chirality structure factors for K — 0.1 (left), 
K = 0.4 (center), and K = 0.7 (right). Inset of the bottom 
panel: A Lorenzian fit of the peak at = n and a; ~ with 
77 = 0.1. The dashed lines denote the lower and upper edges 
of the continuum. 

Qx = and TT. The existence of the gap implies that the 
scalar-chirality order still belongs to an excited state. If 
we assume a complete staggered dimer order, i.e., the 



ground state is the product of local dimcr singlets, the 
scalar-chirality operator ([6]) may be effectively reduced as 

Sx,i-{Sx+is'><Sa;+i,2)\i'o) ~ {Sx,i/2)\'ipo). Thus, the dis- 
persions arc similar to those of the spin structure factor in 
the onc-dimensional spin-Peierls Heisenberg modeli^i^ 
For K ~ 0.7, we can see the closing of the gap and, more- 
over, the appearance of a dominant peak at 5^ = tt and 
Lu ^ 0. This peak is well fitted by a Lorenzian with 
T] = 0.1, as shown in the inset of the right panel of 
Fig. mb). Hence, the spectrum for = tt is composed 
of a (5-function peak at w = and a gapfuU continuum, 
as is D{q,u!) in the staggered dimer LRO phase. It must 
indicate the presence of the scalar-chirality LRO. 

In summary, we study the two-leg Heisenberg ladder 
with the cyclic four-spin interaction. The static and dy- 
namical structure factors for the spin waves, staggered 
dimer order, and the scalar-chirality order parameters are 
calculated with the DDMRG method. We find that the 
spin-spin correlation with wave vector (qx-, qy) = 0) is 
enhanced by the cyclic interaction. We also confirm the 
presence of long-range order in the staggered dimer and 
scalar-chirality phases by a J-function peak contribution 
of the structure factor at a; = 0. 
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